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Abstract 



The electric dipole moment of the muon (cL) is evaluated in a two Higgs 
doublet model with a softly broken discrete symmetry. The leading con- 
tributions from both one loop and two loop diagrams are considered. For 
tan/3 = 1^2 1 /| ^i | close to one, contributions from two loop diagrams involving 
the t quark and the W boson dominate; while for tan /3 ^ 10, contributions 
from two loop diagrams involving the b quark and the r lepton are dominant. 
In these two regions, d^ « (m^/m e )d e . For 8 > tan/3 > 4, significant cancel- 
lation occurs among the contributions from two loop diagrams and the one 
loop contribution dominates for tan j3 ~ 7. For tan f3 > 15, the calculated 
can be close to the reach of a recently proposed experiment at the Brookhaven 
National Laboratory. 
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I. INTRODUCTION 



A non-self-conjugate particle with a non-zero spin can have a permanent electric dipole 
moment (EDM) if the combined transformation of charge conjugation (C) and parity (P) 
is not an exact symmetry. In the Standard Model (SM) of electroweak interactions, CP 
violation is generated by the Kobayashi-Maskawa (KM) phase. The electron EDM generated 
from the KM mechanism is about 8 x 10~ 41 e • cm [|T]-§[. Similarly, the expected muon EDM 
in the SM is about 2 x 10~ 38 e • cm. These values are more than 13 orders of magnitude 
below the current experimental limits. Therefore, precise measurements of the electron and 
the muon EDMs might reveal new sources of CP violation beyond the SM. 

Several experiments have been carried out to search for a neutron EDM (d n ) and an 
electron EDM (d e ). At present, the experimental upper limits on d n [|J and d e [[J at 95% 
C.L. are \d n \ < 11 x 10~ 26 e • cm and \d e \ < 6.2 x 10~ 27 e ■ cm respectively. The measurement 
of the muon EDM (d^) is not yet so precise as that for d e and d n . The experimental upper 
limit on the muon EDM at 95% C.L. is \d^\ < 1.1 x lCT 18 e • cm ||. Recently, a dedicated 
experiment has been proposed Jj]] to measure the electric dipole moment of the muon with the 
existing g — 2 ring at the Brookhaven National Laboratory (BNL) ||. This experiment may 
be able to improve the measurement of the muon EDM by at least 4 orders of magnitude. 

A significant electric dipole moment for the electron or the neutron can be generated if 
CP violation is mediated by Higgs-boson exchange [PHTT||. In a non-supersymmetric model 
with multi-Higgs doublets, for flavor symmetry to be conserved naturally to a good degree, 
a discrete symmetry is usually required fl2[]. In a two Higgs doublet model, there will be 
no CP violation from the Higgs sector if the discrete symmetry enforcing natural flavor 
conservation is exact. If this symmetry is broken only by soft terms, CP violation can be 
introduced and the flavor changing interaction can be kept at an acceptably low level as well 

MM- 



The EDM of an elementary fermion, generated from Higgs-boson exchange, has dominant 
contributions from one-loop or two-loop diagrams. Intuitively, the one-loop contributions 
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are proportional to the product of two Yukawa couplings (m/v) 2 ; while the leading two-loop 
contributions are proportional to only one Yukawa coupling m/v, where m is the fermion 
mass and v = 246 GeV is the vacuum expectation value (VEV) of the SM Higgs field. 
The exact one loop result actually goes like m 3 /m^hi(m 2 /mg), where mo is the mass of the 
lightest neutral Higgs boson. The two-loop contributions dominate for the electron EDM 
because its one-loop contributions are greatly suppressed. For the muon EDM, the one-loop 
contributions could be comparable to those from two-loop diagrams for large tan (5. 

Several studies have confirmed that a significant electron EDM can be generated if Higgs- 
boson exchange mediates CP violation in a two Higgs doublet model with Yukawa interac- 



tions of Model II [|T^,|T6|]. In this Model, one Higgs doublet couples to down- type quarks and 
charged leptons while another doublet couples to up-type quarks and neutrinos. The elec- 
tron EDM has important contributions from two-loop diagrams involving the top quark (11 



and the W boson ||T7Hr9|]. Diagrams with the charged Higgs boson ( J ff ± ) can also contribute 



to the electron EDM . However, recent measurements on the decay rate of b — > by the 
CLEO collaboration |?TJ constrain the mass of the H in the Model II to be at least three 



times larger than the W boson mass and contributions from the H are therefore 

not expected to be significant. For the muon EDM, the two loop contributions involving 
the t quark as well as the W boson were found to be dominant and to satisfy a simple 
scaling relation d^/de ~ m IM /m e [p8,2"9f]. However, this conclusion did not take into account 



potentially important two loop contributions involving the 6-quark and the r-lepton. 

In this paper, we present the first complete calculation for the muon EDM in a two 
Higgs doublet model with CP violation generated from Higgs-boson exchange and Yukawa 
couplings of Model II. We evaluate the leading contributions to the EDM for the muon from 
two loop diagrams involving the W boson, the t quark, the h quark and the r lepton and the 
one loop diagrams. We find that (i) for tan/5 > 10 contributions from two loop diagrams 
involving the h quark and the r lepton are dominant ||30|| ; (ii) for 8 ^ tan/3 ^ 4, significant 
cancellation occurs among the two loop diagrams and the one loop contribution dominates 
for tan (3 ~ 7. 



II. CP VIOLATION FROM HIGGS EXCHANGE 



In a two Higgs doublet model with a discrete symmetry softly broken, the Higgs potential 



can have the following form [10| 



V[(f>l, 2 ] = ^10101 + "120202 + #102 + ??*0201 
+ 7^l(0 t 10l) 2 + ^2(0202) 2 

+S(0 t 10l)(02> 2 ) + ^'(0 t 102)(020l) + \kA<\>2? + ^*(02>l) 2 . 
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This potential respects a discrete symmetry X — > — 0i, and 2 — > +02, except for the soft 
terms 0102 and 0102, We can rewrite the Higgs potential as 
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where all the coupling constants are real, with the possible exception of h. CP is violated if 
hv* 2 v\ has an imaginary part. The minimum of this potential occurs at 



< 01 >= 



y/2' 



< 02 > = 



V2 
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where v\j y/2 and V2/V2 are the vacuum expectation values (VEVs) of 0i and 02- Both VEVs 
can be complex. Without loss of generality, we will take < 0^ >= Vi/y/2 and < 0° >= 



t^e* 6 * I y/2 with V\ and v-i real and tan/5 = f 2 /Vi- The two VEVs satisfy yvf + v^ — v, where 
v is the VEV of the SM Higgs field and v 2 = (V2G F )-\ 

Introducing a transformation, which takes the two Higgs doublets to the eigenstates $1 
and $ 2 , such that < $5 >= v/y/2 and < $^ >= 0, we have 



!>i = cos/9$i - sin/3$ 2 , 02 = (sin/?$i + cos/3$ 2 )e 
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where G ± and G° are Goldstone bosons, H ± are singly charged Higgs bosons, H\ and H 2 
are CP-even scalars, and A is a CP-odd pseudoscalar. 

In the new eigenstates of Eq. [| for the Higgs fields, the Higgs potential becomes f20|| , 
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$ 2 ] = ^Al($i$! - y ) 2 + -A 2 ($5$ 2 ) 2 

+ A 3 ($I$i - y)$£$2 + A 4 ($ t 1 $ 2 )($^ly 
V 2 

+ A 5 ($ t 1 $i + $£$ 2 - y)($I$2 + $2$l) 



+ (A 6 $ t 1 $ 2 + A*$5$ 1 )($ T 1 <I> 1 - $3,$ 2 - y) 

+^A 7 ($ t 1 $ 2 ) 2 + iA;(^$ 1 ) 2 
+p($^$ 2 ), 
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where the parameters p, t> and Ai, % = 1 through 5, are all real; Xq and A7 can be complex. 
CP is violated if the imaginary part of Xq or A7 is nonvanishing. There are 11 parameters, 
but only 10 of them are independent. One of the relations among the 11 parameters is 

As 



Im(A6) = -2(A 1 -A 2 ) 
In the above parameterization, tan j3 is given by 

Ai — A2 
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In Model II [|15|,|T6|], the down- type quarks and charged leptons have bilinear couplings 
with the doublet <f)\, while all the up-type quarks and neutrinos have bilinear couplings with 

the doublet 2 . The Lagrangian density of Yukawa interactions has the following form 

3 3 
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l m , d m , and w m are the leptons, the down-type quarks and the up-type quarks in the gauge 
eigenstates. This Lagrangian respects a discrete symmetry, 



Lm 
L 



' b Ri a R ~ * a Rl 



(ii) 



with m = 1, 2, 3 and a = 1, 2. 

In the new eigenstates of Eq. [| for the Higgs fields, the neutral Yukawa interactions of 
the quarks become 



J2 —ll(H 1 -taa/3H 2 )-i ^ — h 5 l(G° - tan (3 A) 
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u=u,c,t 



where the quarks and leptons are in the mass eigenstates. 

Adopting Weinberg's parameterization [IU| we can write the following neutral Higgs- 
boson exchange propagators as 

-,2 a v i ™„2 
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where the summation is over all the mass eigenstates of neutral Higgs bosons. We approxi- 
mate the above expressions by assuming that the sums are dominated by the lightest neutral 
Higgs boson of mass mo, and drop the sums and indices n in Eq. [L3] hereafter. There are 
relations among the CP violation parameters: 



ImZo + ImZo = — cot 2 /JImZi, 
lmZ — lmZ = + tan 2 /3ImZ 2 . 



(14) 



Employing unitarity, Weinberg has shown [|1(J that |ImZi| < (l/2)tan/3(l + tan 2 /5) 1//2 and 
\ImZ 2 \ < (l/2)cot/3(l + cot 2 /3) 1/2 . Therefore, we have 



\lmZ + lmZ \ < (l/2)cot/3(l + tan 2 /3) 1/2 , 
\lmZ -lmZ \ < (l/2)tan/3(l + cot 2 /3) 1/2 . 



(15) 



III. MUON ELECTRIC DIOPLE MOMENT 



A. Two Loop Diagrams 



Two loop diagrams for fermion loops contributing to the muon EDM are illustrated in 
Figures 1(c) and 1(d). The diagrams with the intermediate Z boson are highly suppressed by 
the vector part of the Zp + p~ couplings. Therefore, we consider only the diagrams involving 
an intermediate 7. 



The muon EDM generated from two loop diagrams with the top quark is [11 



ce loop = -™ emaV * GF {[f( Pt ) + g(jh)]1mZ + [g( Pt ) - f( Pt )]lmZ Q }, 



(16) 



3 (4vr) 3 

where m = m^, p t = mf/rriQ, and mo is the mass of the lightest neutral Higgs boson. The 



functions / and g are defined as 
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We take m t = 175 GeV, m b = 4.8 GeV, m T = 1.777 GeV, m w = 80.0 GeV, and the fine 
structure constant a = 1/137, neglecting the running up to the scale of mo- 
The EDM generated from the b and the r loops is [|30| 
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^r-ioop = _ (4 iV e g 2 tan 2 p) Cma ^ F {f(p f ) + g(PfWmZ + ImZ ), 



(18) 



where N c is the color factor and Q is the charge. For the b and the r, 4N C Q 2 is equal to 4/3 
and 4 respectively. 



The leading contribution from two loop diagrams with the W boson fl7j-|19 
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(19) 
(20) 
(21) 

x(l — x) — r ( r J ' ^ ^ 

where p^ = ^w// m o> an d the functions / and g are defined in Eq. [T7]. There are another two 
sets of diagrams with the W boson [T7[ contributing to the muon EDM. The contribution 
from these additional diagrams has an opposite sign to that of dY~ l °° v and therefore reduces 
the magnitude of the full contribution from the W boson. In our analysis, we have employed 



r f 1 1 
h(r) = - / dx— 



x) — r 



the formulas in Ref. 17 to evaluate the complete contributions from the W loops. 



B. One Loop Diagrams 

The one loop diagrams contributing to the muon EDM are illustrated in Figures 1(a) 

and 1(b). The one- loop contribution to the muon EDM (o^) is 

,i-ioop emV2G F tan 2 P f ~ 

p = ^y 2 I(p)(ImZ + lmZ ), (23) 

where p = m? /m^. The function I(p) is defined as 

i x 2 



I(p) = p I dx 



o px 2 — x + 1 

= 1 + f 1 J (*i " 1) ) - (to - 1) M^)], (24) 

p{z\ - Z 2 ) Zi z 2 

where z\ and z 2 are roots of px 2 — x + 1 = 0. For p « 1, I(p) approaches — p[ln(p) + 3/2]. 
In this limit, the one-loop contribution to the the muon or the electron EDM has a simple 
form 
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d\$* = - 6V (/^ n ^ m 3 /m 2 )[Hm 2 /ml) + 3/2](ImZ + lmZ ). (25) 

where m = me^m^ and (mo >> m) is assumed. Therefore, the one-loop contribution to the 
muon or the electron EDM is proportional to (m 3 /mg) ln(m 2 /wig). 



C. Numerical Values 

In this section, we discuss the numerical value of the muon EDM with several choices of 
mo and tan/3. The contributions from one loop diagrams (1-loop) and two loop diagrams 
involving the W boson (W-loop), the t quark (t-loop), the b quark (6-loop) and the r lepton 
(r-loop), are presented in Table I for tan/3 = 2 and 20 with m = 100, 200 and 400 GeV, in 
units of (a) Im Zq and (b) Im Zq. 

To study the muon EDM dependence on the lightest Higgs boson mass (mo), we present 
the muon EDM generated from one loop and two loop diagrams in units of (a) Im Zq and 
(b) Im Zq, as a function of mo, with tan/? = 1, 7 and 20 in Figs. 2, 3 and 4. Figure 5 shows 
the effect of varying tan/? on the muon EDM, for the case of m = 100 GeV. From these 
Figures, we conclude that (i) for tan/3 close to one, contributions from two loop diagrams 
involving the top quark and the W boson dominate; (ii) while for tan/3 > 10 contributions 
from two loop diagrams involving the b quark and the r lepton are dominant ]30|]; (hi) 
for 8 ^ tan/3 ^ 4, significant cancellation occurs among the contributions from two loop 
diagrams and one loop diagrams dominate for tan ~ 7. 

To compare the muon EDM with the electron EDM, we present the ratio of (d^/de) to 
(m^/m e ) in Fig. 6. For the electron EDM, the two-loop contribution is about 10 5 (for tan (3 > 
10) to 10 6 (fortan/? ~ 1) larger than that from the one-loop diagrams. For the muon EDM, 
the one loop diagrams make important contributions which dominate for tan/3 ~ 7. For 
tan/3 > 10, the one-loop contribution to the muon EDM is about 10% of that from the 
two-loop diagrams. We find that (i) for tan/3 ~ 1, d M « (m^/me); (ii) for tan/3 ^ 10, 
dp ~ 0.9(m^/m e ), since there is a cancellation for the muon EDM between the one-loop 
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contribution and the two-loop contribution involving the 6-quark and the r-lepton; (iii) for 
8 > tan/5 > 4, \d^\ can be two to three times \(m^/m e )d e \. 
The muon EDM can be expressed as 

dfj, = d A lmZ + d B lmZ 

= ^±^(lmZ + lmZ ) + ^^(lmZ - lmZ ) (26) 

where d A and d% are the total coefficients of lmZ and lmZ from one loop diagrams, the 
H^-loop, the i-loop, the 6-loop, and the r-loop. Applying unitarity constraints in Eq. [15|, we 
can define the maximal \d^\ as 

Wmax = l^^lilmZo + ImZol + |^%mZ Q - ImZ | 

= \^A±M cot + tan 2 (3) 1 ' 2 + t an/5(l + cot 2 (3) 1 ' 2 . (27) 

In Figure 7, we present the maximal value allowed by unitarity for the muon and the 
electron EDMs, as a function of tan/5 with m = 100, 200 and 400 GeV. Also shown is 
the experimental upper limit for the electron EDM (denoted by |rf e |ux.)- I n this model, for 
tan/5 > 10, the simple scaling {jn^/m e )\d e \-[j^ could be treated as an upper limit for the 
muon EDM since d^ ~ (m M /m e )c? e . A measured value of the muon EDM above this bound 
could arise if the muon EDM and the electron EDM are generated from different sources, i. 
e., if this model applies only to the muon EDM. 

There are several interesting aspects to note from the different contributions: (i) The 
contributions from one loop diagrams and two loop diagrams with the b and the r are 
proportional to (Im Z +Im Z ). (ii) For the the t-loop, the coefficient of the Im Z is much 
smaller than that of the Im Zq. (iii) The VT-loop does not contribute to the Im Zq term, 
(iv) The contributions from one-loop diagrams and the PF-loop have positive sign while the 
contributions from heavy fermion loops have negative sign. Therefore, the total muon EDM 
is slightly reduced by cancellation, (v) For large tan/5, two loop diagrams involving the b 
and the r dominate and the electron EDM and the muon EDM are almost proportional to 
their masses. Therefore, d M ss {m^/m e )d e for tan/5 > 10. 
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IV. CONCLUSIONS 



Our results may be summarized as follows: 

• For tan /3 close to one, contributions from two loop diagrams involving the t quark and 
the W boson dominate. 

• For tan j3 > 10, contributions from two loop diagrams involving the b quark and the r 
lepton are dominant. 

• For tan/3 ~ 1 or tan/3 > 10, ~ {m^/m e )d e . 

• For 8 ^ tan (3 ^ 4, significant cancellations occur among the contributions from two 
loop diagrams and the one loop contribution dominates for tan (5 ~ 7. 

• For tan/3 > 15 and the lightest Higgs boson mass mo < 300 GeV, CP violation 
mediated by Higgs-boson exchange in a two Higgs doublet model could produce a 
muon EDM which is close to the reach of the proposed BNL experiment 0. 

High precision experimental measurements of the muon EDM could provide to inter- 
esting information about mo and tan/3 as well as CP violation parameters, ImZo, and 
ImZi, i = 0,1,2. A positive result for the muon EDM measurement could shed light on non- 
standard CP violation and help pin down the value for tan /3 in two Higgs doublet models. A 
negative result, however, could mean (i) that tan/3 is small than 10; or, (ii) the CP violation 
parameters ImZ , and ImZj are smaller than their unitarity bounds; or, (iii) the masses of 
the neutral Higgs scalars and the Higgs pseudoscalar are very close to one another (K| . 
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TABLES 

TABLE I. The muon EDM from one-loop diagrams (1-loop) and two-loop diagrams with the 
W boson (PL-loop), the top quark (t-loop), the bottom quark (6-loop), the r-lepton (r-loop) and 
their total (Total) for tan [3 = 2 (tan (3 = 20) with the lightest Higgs boson mass mo = 100, 200 
and 400 GeV, in units of (a) ImZ x 10~ 24 e • cm and (b) ImZ x 10~ 24 e • cm. 







m (GeV) 




Diagrams 


100 


200 


400 


(a) Im Z x 10" 24 e • cm 








1-loop 


+0.01 (+1.19) 


+0.003 (+0.33) 


+0.001 (+0.09) 


TT-loop 


+2.47 (+3.07) 


+1.16 (+1.45) 


+0.20 (+0.25) 


i-loop 


-1.88 (-1.88) 


-1.23 (-1.23) 


-0.72 (-0.72) 


6-loop 


-0.06 (-5.61) 


-0.021 (-2.08) 


-0.007 (-0.73) 


r-loop 


-0.04 (-3.97) 


-0.014 (-1.36) 


-0.004 (-0.44) 


Total 


+0.50 (-7.2) 


-0.10 (-2.8) 


-0.53 (-1.5) 


(b) Im Z x 10~ 24 e • cm 








1-loop 


+0.01 (+1.19) 


+0.003 (+0.33) 


+0.001 (+0.09) 


VL-loop 


+0.0 


+0.0 


+0.0 


i-loop 


-0.33 (-0.33) 


-0.21 (-0.21) 


-0.11 (-0.11) 


6-loop 


-0.06 (-5.61) 


-0.021 (-2.08) 


-0.007 (-0.73) 


r-loop 


-0.04 (-3.97) 


-0.014 (-1.36) 


-0.004 (-0.44) 


Total 


-0.42 (-8.7) 


-0.24 (-3.3) 


-0.13 (-1.2) 
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FIGURES 

FIG. 1. Feynman diagrams for leading contributions to the muon EDM from one loop diagrams 
[(a) and (b)] and two loop diagrams with heavy fermions and the W boson [(c) and (d)]. There 



are many more diagrams involving the W boson |17-|l9f| that are not shown in this figure. 



FIG. 2. The muon EDM generated from one-loop diagrams (dash), two-loop diagrams 
(dash-dot) and their total (solid) in units of (a) Im Zq and (b) Im Zq, as a function of mo, 
with tan/3 = 1, where mo is the mass of the lightest neutral Higgs boson. Also shown are the 
contributions from two- loop diagrams involving the W boson (dot), the t quark (dash-dot), the b 
quark (dash-dot-dot), and the r lepton (dash-dot-dot-dot), in units of (c) Im Zq and (d) Im Zq. 



FIG. 3. The same as in Fig. 2, except that tan/3 = 7. 



FIG. 4. The same as in Fig. 4, except that tan/3 = 20. 

FIG. 5. The muon EDM generated from one-loop diagrams (dash), two-loop diagrams 
(dash-dot) and their total (solid) in units of (a) Im Zq and (b) Im Zq, as a function of tan/3, with 
itiq = 100 GeV. Also shown are the contributions from two-loop diagrams involving the W boson 
(dot), the t quark (dash-dot), the b quark (dash-dot-dot), and the r lepton (dash-dot-dot-dot), in 
units of (c) Im Zq and (d) Im Zq. 

FIG. 6. The ratio of (d^/de) to (m M /m e ) in units of (a) Im Zq and (b) Im Zq, as a function 
of tan/?, for m = 100 GeV (solid), 200 GeV (dash) and 300 GeV (dash-dot). 

FIG. 7. The maximal EDM allowed by unitarity (Eq. [27] ) for (a) the muon EDM and (b) 
electron EDM, as a function of tan/3 for mo = 100 (dash), 200 (solid) and 400 GeV (dash-dot). 
Also shown is the experimental upper limit (U.L.) on the electron EDM at 95% C.L. (dot). 
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